Pairing of a trapped resonantly-interacting fermion mixture with unequal spin 

populations 
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We consider the phase separation of a trapped atomic mixture of fermions with unequal spin 
populations near a Feshbach resonance. In particular, we determine the density profile of the two 
spin states and compare with the recent experiments of Partridge et al. Q. Overall we find quite 
good agreement. We identify the remaining discrepancies and pose them as open problems. 



Introduction — The usual Bardeen-Cooper-SchriefTer 
(BCS) theory of Bose-Einstein condensation (BEC) of 
fermion pairs requires the populations of the two species 
involved in the s-wave pairing to be equal. For a long 
time, therefore, theorists have discussed fermionic pair- 
ing when the species densities are unequal, and several 

Soposals for the ground state have been put forward 
Experimentally, however, such superfluid states 
with unequal densities have remained elusive. 

After several years of experimental studies of the BEC- 
BCS crossover with equal spin population, experiments 
with ultracold atoms have very recently also turned to 
studying superfluidity with unequal populations 0, 
The basic idea is to load a trap with an unequal pop- 
ulation of two hyperfine states of ^Li and tune the bias 
magnetic field close to a Feshbach resonance. It turns out 
that the physics of an unequal-population Fermi mixture 
in a trap is rather different from the uniform case, the 
dominant characteristic of the measured density profiles 
appears to be a phase separation, with an equal-density 
BCS phase in an interior core, and an outer shell consist- 
ing mostly of the majority species. If there are additional 
two-component phases, for instance Fulde-Ferrel-Larkin- 
Ovchinikov (FFLO) phases, they are confined to a shell- 
shaped region between the outer majority shell and the 
inner BCS-like core. 

Partridge et al. have reported in-situ measurements 
of the density profiles of the two states Q, while Zwier- 
lein et al. report on measurements after expansion [^. 
The former experiments are performed close enough to 
the Feshbach resonance that they may be regarded as 
being in the unitarity limit, i.e., in the limit that the in- 
teraction strength g is effectively infinite so it does not 
provide an energy scale to the problem. In this Letter we 
concentrate on the data from Ref. and limit ourselves 
therefore to the unitarity region. We do not consider here 
the data of Ref. j3] on rotating fermion gases, nor do we 
deal with the issues arising from expansion after the trap 
is switched off. 

We present a zero-temperature analysis of the phase 
separation using a local density approximation (LDA) 
and a BCS ansatz for the many-body wavefunction. The 
Feshbach resonance is treated using a single-channel de- 



scription, because the closed-channel component of the 
Cooper pair wavefunctions is small in the crossover re- 
gion for the extremely broad Feshbach resonance that is 
being used in the experiment Q . Based on an anal- 
ysis of the uniform case at unitarity, we give simple ar- 
guments for the occurrence of phase separation, and to 
identify the surface that surrounds the BCS phase. We 
then calculate the majority and minority density profiles 
within the BCS ansatz, and compare with experimental 
profiles. 

BCS ansatz for the unitarity regime — We first ex- 
amine pairing at unitarity with unequal chemical poten- 
tials for a homogeneous mixture. Since we will treat the 
trapped case in LDA, the results from this analysis can be 
used locally for any point in the trap. We are interested 
in the g —t oo limit of the Hamiltonian 



H 
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The index a runs over the two hyperfine states of ^Li, 
denoted by |1) and |2). The masses are the same, so ek = 
fi^V^ /2m for both species, but the chemical potentials are 
different, i.e., /ii^2 = /i ± ft, , so that it is possible to have 
unequal densities ni^2 = n ±m . 

We use the BCS wavefunction as an ansatz for 
the paired ground state. This corresponds to using 
the following decomposition for the interaction term: 



''k,l''-k,2 



+ C_k,2Ck,l 



A'^/g. We restrict our- 



selves to zero-momentum pairing, because the experi- 
mental data do not indicate the presence of an FFLO 
state and because two-channel calculations suggest 
that FFLO states are not stable close to the resonance 
At unitarity, the BCS ansatz is best understood as a 
variational approach as opposed to a mean-field approx- 
imation. For the case of an equal-density mixture, this 
approach has been shown to be a reasonable method of 
interpolating between the BCS and BEC limits. We em- 
ploy the same philosophy here and extract information 
using the BCS expressions for the number density, energy 
density and the gap equation. Fortunately, improved in- 
formation about the equal-density case is available from 
Monte-Carlo simulations [liillllllJI and may be used to 
improve our trap calculations. 
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For the equal-density case, and within the single- 
channel assumption, the many-body system at resonance 
is universal in the sense that the only energy scale in the 
problem is that set by the density, i.e., the Fermi energy 
of the corresponding free gas ep- The BCS ansatz gives 
the energy of the resonantly interacting system to be 0.59 
times ep . This number is called 1 + f3, with the universal 
number (3 being known from Monte Carlo calculations 
to be ^ ~ -0.58 US [H E3. to which the BCS value 
/3 ~ —0.41 should be regarded as an approximation. In 
the strong-coupling region, the pairing gap is of the order 
of the chemical potential, instead of being exponentially 
suppressed as in the weak-coupling regime. Within the 
BCS ansatz Aq — 1.16/z ~ 0.68eF, while Monte-Carlo 
calculations give Aq ~ 0.84eF [l^ . 

With differing chemical potentials, i.e., h ^ 0, the 
quasiparticle energy spectrum of the BCS ansatz has two 
branches Ei,^±^ Ek±h, with = ^/A^T^ 0. For 
h > A, the lower branch becomes negative in the mo- 
mentum interval (fci, fc2), with ki,2 = /i =F v/i^ — A^. In 
that case we need to fill up the negative energy modes 
to construct the lowest-energy ground state. As a result, 
the thermodynamic potential becomes 



V 4^ 



(fk - ^J■- Ek) 
+ eih-A) 
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fci<|k|<fe2 



(1) 



The A^/2ek and A'^ /g terms are required to remove the 
usual ultraviolet divergence. Note that the last term is 
the only place where h enters. The gap equation is given 
by the condition for the extrema of the thermodynamic 
potential &{fi,h,A) i.e., dQ/dA = 0. The density av- 
erage and difference may be calculated as n = —d0/dfj, 
and m = —d0/dh. The resulting expressions also have 
contributions from the (fci,fc2) shell in addition to the 
usual BCS contributions. 

The ground state of the system is the absolute mini- 
mum of 0. For equal chemical potentials, the function 
0(/z,O, A) has a maximum at A = and a minimum 
at the equal-density gap Aq. As h is increased, there 
is a certain value, h = Ao/ai, at which the A = ex- 
tremum becomes a minimum and there is an intermediate 
maximum. This maximum corresponds to an unequal- 
density paired solution of the gap equation, commonly 
known as the Sarma phase 0, which is thus unstable. 
At some higher value of the chemical potential difference, 
h = Ao/a2, the A = solution becomes the global min- 
imum so that the normal state is more stable than the 
paired state. For weak coupling {g^^ — oo) the special 
values of /i are h = A/2 and h = A/\/2 ~ A/1.414. Mov- 
ing towards the resonance from the BCS regime, we find 
that the values of ai and a2 change only slightly within 
the BCS ansatz. The change in ai is smaller than 1%, 
while we find that a2 evolves to about 1.44 at unitarity, 
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FIG. 1: Left: argument for phase separation from consid- 
eration of the local effective chemical potentials. Shown are 
fJ'csir) (dashed) and fics{r) ± h (solid). Right: phase separa- 
tion, after the x and y coordinates have been scaled to make 
the trap look spherically symmetric. Here A is the BCS core, 
C is the outer majority shell, B is the intermediate shell of bi- 
component phase which we treat as two non-interacting ideal 
gases. 



i.e., a change of about 2%. The value of a2 at infinite 
coupling is a universal number and our value a2 ~ 1.44 
is the approximation to this universal number within the 
BCS ansatz. 

Knowing that the Sarma phase is a maximum of the 
thermodynamic potential in the uniform case, it is clear 
that we will not have this phase in the trap within the 
local-density approximation. Thus, barring more exotic 
pairing mechanisms, we only have to consider normal 
phases and equal-density BCS phases. 

Phase separation in a trap — The trapping poten- 
tial used in Ref. j^] is asymmetric and obeys 14rap(r) = 
^m0lz'^ + \m0\{x'^ + y'^) , with 0^ = 2'kx (7.2 Hz) and 
0j_ = 27rx (350 Hz). We scale the spatial variables in 
the radial directions so that the trap potential becomes 
spatially symmetric, with trapping frequency 0z in each 
direction. In LDA, the trap terms in the Hamiltonian are 
absorbed into the chemical potential, so that we have ef- 
fective space-dependent chemical potentials: 



IJ-isi"^) = Ail - Hrap(7-) = (/-i + h)- \m0lr 



(1') 

Moff (^) = - Vtrap(f-) ^(p~h) 



w2 2 

m0,r 



(2) 



The average ^^ff — ^ ^Moff + t^csj decreases parabol- 
ically away from the center of the trap while the dif- 
ference equals 2h and stays constant. Near the cen- 
ter of the trap, h is small compared to ficS and hence 
compared to Aq{^cs)- Thus the densities are forced 
to be equal and we have a BCS phase in the center. 
Since /icff(r) decreases monotonically, there is some ra- 
dius i?BCS at which Ao(/ioff) — 1.16/ioff is equal to a2h. 
Outside this radius, a two-component normal phase is 
more stable than a superfluid state. For this phase, 
we ignore interactions between the two components and 
treat it like two ideal Fermi gases whose densities are 
determined by their different chemical potentials. At 
the Thomas- Fermi radius of the minority (|2)) species, 
R^p = \/2{p — h)/m0'^, the minority density vanishes. 
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FIG. 2: Majority and minority densities calculated in LDA 
for typical experimental parameters of Ref. 3- The densities 
are expressed in units for the effective spherical trap with 
trapping frequency 0z in each direction after rescaling the 
radial directions. 

(2) 

Outside r — Wj^p only the majority remains. This outer 
shell survives up to the majority (|1)) Thomas- Fermi ra- 
dius R^r^l = v/2(^ + h)/m0l. 

Density profiles — Given the chemical potentials fi±h 
at the trap center, we can calculate the BCS core radius 

..(12) 

_Rbcs a-nd the Thomas-Fermi radii Rj^p . The densities 
are then given by 

(nBcsilJ-csir)) r < Rbcs 
mifJ-cff{r) ± h) i?BCS <r < i?^^/^ (3) 
r> i?^^/^ 

Here ^bcsIa^) = Z^k [1 ^ ('^k - fJ-)/Eu] /2V is the usual 
equal-density BCS density for a single component, and 
?in(m) = (2m/i)^/^/67r^ is the ideal-gas density. At uni- 
tarity, n-BcsifJ-) = "n(m/[1 +/?])■ To calculate the density 
corresponding to given total atom numbers, we first find 
the chemical potentials which give Ni^2 = J drni^2{'''), 
and then use the expressions given above. 

In Fig. 12 we show density profiles for a typical experi- 
ment of Ref. . Experimentally, the density itself is not 
accessible, but it is interesting to note some features. At 
the BCS core radius -RbcS: the majority density ni{r) 
has a discontinuity because it is determined by the BCS 
density corresponding to ^ on the r < Rbcs side and 
by the normal density corresponding to n + h on the 
r > i?BCS side. In the weak-coupling limit, the nBcs(M) 
and TIN (a*) functions are almost identical, and so the 
discontinuity would then have been much more promi- 
nent. At unitarity, however, nBCs(M) = "■n(/^/[1 + P]) = 
[1 + l3]^'^/^n]<s{^), so that n]<s{^) ~ 0.454riBcs(A^) for our 
BCS treatment. This significantly reduces the upward 
jump of the majority density at the core edge. 

Similarly, the minority density n2 (r) has a discontinu- 
ity at the BCS core edge as it jumps down from nBcsifJ-) 
to nT<i(fj,— h). This discontinuity is enhanced by the effect 
of nonzero (3 in n]<s{^) = [l+/3]^^^nBcs(M)- The large de- 
crease of the minority density assures that the minority 

(2) 

Thomas- Fermi radius i?rpp is only slightly larger than the 
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P = (N -N,) / (Nj+N^) 
FIG. 3: Radius versus number asymmetry. The three solid 

(2) 

lines from bottom to top are calculated values of iiscs, i?TP 
and R^Yp. Filled circles and empty squares are experimental 
majority and minority radii from Ref. Q. The radii Rbcs, 

(2) 

7?^p and the experimental minority radii are scaled by the 
ideal-gas Thomas- Fermi radius R^p^^(N2) of N2 (minority) 
fermions. The radius rI^^ and the majority radii are scaled 
by 



BCS core radius i?BCSi i-e., that the intermediate two- 
component shell is rather thin. In the real system, we ex- 
pect the LDA discontinuities to be smoothed out some- 
what by gradient and other corrections. Since experi- 
mentally only spatially integrated versions of ni(r) (col- 
umn densities) are observed, the small non-monotonicity 
of the majority density is further washed out and is ex- 
pected to be difficult to observe. 

Majority and minority radii — Fig. 13 shows the evo- 
lution of the three radii in our theory (i?BCSi -^tF' ^tf) 
with the number asymmetry P — {Ni — N2)/{Ni + N2), 
and compares with radii from Ref. jj]. Measured in 
units of the ideal-gas Thomas-Fermi radii correspond- 
ing to Ni^2, the theoretical curves depend only on P and 
not on the total number. The experimental radii are ex- 
tracted by fitting the measured column density profiles 
to ideal-gas Thomas-Fermi distributions. 

It is reasonable to assume that the experimental minor- 

(2) 

ity radii correspond to i?BCS rather than i?rpp, since the 
minority occupancy in the intermediate shell is negligible, 
as shown in Fig. |21 Our ab initio LDA calculations then 
explain the radius data extremely well at large P. At 
small P, the calculated radii are somewhat higher than 
the experimental ones. This is expected from the use 
of the BCS ansatz, which underestimates the reduction 
of the paired state energy, and hence also the reduction 
of the size in a trap. We could improve our calcula- 
tion by using the Monte Carlo values /.t ~ 0.42eF and 
Ao — 0.84eF. However, to identify i?BCS we also need 
the ratio a2 — Aq/H at which the unequal-density nor- 
mal phase becomes more stable than the equal-density 
paired phase, and to the best of our knowledge 012 is not 
known outside the BCS ansatz. We have therefore opted 
for consistency and used the BCS ansatz throughout. 
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The question remains whether there is a critical 
nonzero value of P at which phase separation first ap- 
pears. No such feature appears in our calculations, be- 
cause with the BCS treatment of unitarity, phase sepa- 
ration appears at any nonzero asymmetry, since only the 
equal-density BCS phase and the normal phase are stable 
in this case. 

Axial density profiles — In Fig. 01 we have plotted 
LDA calculations for densities with both x and y direc- 
tions integrated out, and compared them with column 
densities integrated along the x direction j^. The typi- 
cal feature seems to be that the majority density profile 
fits better than the minority profile, especially outside 
the BCS core. This is not so surprising because, outside 
the core, the majority distribution is simply the ideal-gas 
Thomas-Fermi distribution. Some details of the density 
profiles are smoothed out because of the double integra- 
tion, but it is instructive to look at the axial density dif- 
ference. In the integral for the axial density difference, 
one can for z < Rbcs remove z from both the integrand 
and the integration limits, so that the theoretical axial 
density difference is constant up to z = i?BCS j a-s seen in 
Fig- El In the experimental data, however, there is barely 
an extended constant part in the difference. Geometri- 
cally, the experimental data indicates that the inner core 
is expanded radially and squeezed axially compared to 
the LDA prediction. At this point the reason for this 
discrepancy is not clear. Possible reasons could be tem- 
perature effects, nonuniversal physics beyond the single- 
channel description, or the effects of nontrivial phases 
in the interface region which we have not included. An- 
other intriguing possibility is that, since the trap is much 
tighter in the radial direction than in the axial direction, 
corrections to the local density approximation might be 
required for the radial directions. Chemical potentials 
are typically 5-15 times the energy corresponding to ra- 
dial trapping frequency, the lower end of which might 
be near the limit of validity of LDA. The discrepancy in 
the density difference is an urgent issue and is presently 
under investigation. 

Conclusion — In summary, we have presented ab ini- 
tio calculations of the density profiles for a fermion mix- 
ture near a Feshbach resonance loaded into a trap and 
with unequal spin-populations. While the major features 
are successfully reproduced, two major questions emerge 
from our analysis. One is the shape of the axial density 
difference curve which deviates somewhat from the LDA 
calculation, as seen clearly in Fig. 01 The second issue is 
the possibility of having a transition from a non-phase- 
separated configuration to a phase-separated configura- 
tion at a certain critical value of the number asymmetry 
P. 

We have assumed that unequal-density pairing 
schemes are less favored than the two-component normal 
phase of the intermediate shell. While it is likely that 
the intermediate shell is too small to make a difference 
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FIG. 4: Experimental axial densities as function of z for 
several {Ni,N2), plotted together with theoretical densities 
integrated over both x and y directions: J dxdyn{x, y, z). Up- 
per panels are the axial densities of states |1) and |2) and the 
lower panels are the differences. Note that there are no fitting 
parameters. 



in the density profiles, the issue of stability of various 
unequal-density pairing schemes has not been examined 
thoroughly at unitarity. Some of these questions and is- 
sues are currently also under investigation. 
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Note: At the last stages of our work we learned of 
independent work discussing issues similar to ours [l^ 
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